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Abstract 

We consider classical potential scattering. If at energy E no orbit is 
trapped, the Hamiltonian dynamics defines an integer-valued topological 
degree d.eg(E') < 1. This is calculated explicitly for all potentials, and 
exactly the integers < 1 are shown to occur for suitable potentials. 

The non-trapping condition is restrictive in the sense that for a bounded 
potential it is shown to imply that the boundary of Hill's Region in con- 
figuration space is either empty or homeomorphic to a sphere. 

However, in many situations one can decompose a potential into a sum 
of non-trapping potentials with non-trivial degree and embed symbolic 
dynamics of multi-obstacle scattering. This comprises a large number of 
earlier results, obtained by different authors on multi-obstacle scattering. 

Mathematics Subject Classification: 57R19, 57R20, 70F05, 70F16 



1 Introduction 

In potential scattering on M. d one considers the solutions of the Hamiltonian 
equations for the Hamiltonian function H(p,q) = ^\\p\\ 2 + V(q), where the 
potential V G C 2 (R. d , R) decays at spatial infinity (see (12. 2 ft below), for positive 
values E of H. Equivalently one considers the solutions of Newton's equation 
q = — W(q). The orbit through x is called scattering if rimi t |_ j . 00 \\q(t, x)\\ = oo. 

Here we consider energies E for which there are no trapped orbits, i.e. solu- 
tions where the above limit exists only in one time direction. Then asymptotically 
the solutions have the form of straight lines and can thus be parametrized by a 
point in the cotangent bundle iV := T*S d ~ l . Dynamics induces a diffeomorphism 

S E : N -f N. 



*Mathematisches Institut der Universitat Erlangen-IMiirnberg. Bismarckstr. 1 1/2, D-91054 
Erlangen, Germany, e-mail: knauf@mi.uni-erlangen.de 



1 



In |Kn] this scattering map was used to define a topological index, deg(E) E Z. 
In examples of centrally symmetric V all values < 1 were shown to occur. 

In Sect. |2] we begin by defining the class of potentials for which we can 
explicitly calculate the index in the non-trapping case. This comprises nearly all 
potentials, see Thm. 13. ll for bounded potentials and Thm. I4.3l for potentials with 
singularities, like the Coulomb potential. Billiards can be treated by the same 
method, see Remarks 13.31 

The index is related to the way a Lagrange manifold folds over configuration 
space (Thm. l5~lj) . This shows that only the values < 1 can occur. 

Finally, in Sect. E] we find in all dimensions d topological obstructions for 
(V, E) to lead to motion without trapping. 

In |Kn] this index was used to imbed symbolic dynamics for scattering in a 
potential V = V\ + . . . + T4 for energy E, where the Vi were only assumed 
to carry non-zero degree, and to have non-shadowing supports (no line meeting 
more than two supports). More precisely, for any bi-infinite sequence a in 

{a E {l,...,k} z | aj^ a m }, (1.1) 

there exists an orbit of energy E, visiting the supports of the Vi in the succession 
prescribed by a. So the flow has positive topological entropy if k > 2. 

With the present work, we need not assume any more that the building blocks 
\4 are centrally symmetric, in order to calculate their degree and to combine them 
as indicated above. 

The phenomenon of trapping by chaotic repellers has been observed and 
analyzed in many cases, see e.g. Rapoport and Rom-Kedar [RR] and references 
cited therein. The present work, together with |Kn] , provides a unifying approach 
to several of these results. 

Trapping plays a major role in semiclassical quantum mechanics, and leads to 
the phenomenon of resonances. See e.g. Castella, Jecko and Knauf [CJK], and 
references cited therein. 

Acknowledgement: We thank Christoph Schumacher (Erlangen) and the anony- 
mous referees for helpful comments. 

2 Scattering for Long Range Potentials 

We start by introducing the notions of potential scattering and defining the 
topological degree. 

The configuration space of the scatterer is W 1 , but due to singularities the 
domain of definition M of the potential V may be smaller. 
For physical and mathematical reasons we consider potentials V E C 2 (M, R) for 
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M := R d and M := R d \{s} for some s eR d (in Remarks |Z! we also consider 
the case of several singularities). 

In the 'singular case M = M. d \ {s} we assume that for some Z > 0, a > and 

W G C 2 (M d ,M) 

V{ q ) = ¥ ^ + W{q). (2.1) 

Similar to Derezihski and Gerard |DG] , Sect. 2.7 the force field F := —W of 
V is assumed to meet the long range estimates for multi-indices m G Nq 

/ sup || d m F(q) || r H dr < oo (H < 1) (2.2) 

for some i? (say R = in the non-singular, i? = 2||s|| in the singular case). 

Remarks 2.1 1. For d > 2 the long range condition (12.20 implies the ex- 
istence of limiui-^oo V^g), which we assume to be zero without loss of 
generality. 

Evaluated for m = 0, condition (12. 2D leads to finite total change of velocity 
of the scattered particle (see Thm. 2.5.2 of [DGj). If one would want to 
define so-called M0ller transformations, comparing the dynamics with the 
one for V = 0, a short range condition 

POO 

/ sup \\d m F(q)\\r lml+1 dr < oo (|m| < 1) (2.3) 

JR \\q\\>r 

would be needed (see Thm. 2.6.1 of |DG] ). 

But here we neglect time parametrization of the orbits and consider scat- 
tering on a reduced phase space N. 

2. For all values a > the potential q i— > N j^ |c , meets the long range 
condition ( 12.21) . and for a > 1 the short range condition (12. 3D . 

But for all values a > 2 under the influence of this potential the set of 
initial conditions leading to a collision with the singularity in finite time has 
positive Liouville measure, see, e.g. [LL], §18. We thus assume a G (0, 2). 

We now consider the Hamiltonian function 

HEC 2 (T*M,R) , H(p,q):=l\\p\\ 2 + V(q) (2.4) 

on the symplectic manifold (T*M, u>q), with canonical symplectic form cu := 

Y%=i dqk/\dp k . 
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• For the non-singular case the Hamiltonian flow generated by H on the 
phase space P := T*M is complete (see, e.g. Sect. 2.2 of [DG]). 

• Likewise in the singular case it is known that precisely for a = 2n/(n+ 1), 
n G N the motion can be regularized. For the case of vanishing additional 
potential term W in (12. ip this was treated by McGehee in [MG]. 

For the general case we obtain in Prop. 14.11 below a complete flow on a 
2<i-dimensional symplectic manifold P which (as a set) equals 

P = T*M0(Rx S^ 1 ). 

Physically most important is the case of Coulomb potentials (n = 1). 
In both cases we obtain a flow $ 6 C 1 (M x P,P), also denoted by 

$*:P^P or (p(t,x),q(t,x)) := $*(ar) (t G K), 
restricting to the energy shells T, E := H~ l (E). Moreover 

Knax := sup V(q) E [0, oo). 

Remark 2.2 In potential scattering, the virial identity 



^ hm2 = Jt M*M*)> = 2 ^ E - V ( q W - <?(*)> W(g(t))> (2.5) 



holds true for any trajectory £ i— > (p(t),q(t)) = x), x)) , with energy 
P := iJ(x) (whenever q(t) E M). For E > as a consequence of (12.20 . there 
exists a w'r/a/ radius R V1I = R vir (E) > R, with 

|V(g)|<£/2 and | (q, W (q)) | < E/2 (\\q\\ > R vir ). (2.6) 

Then by O and flU} 

|(g(0,p(0)>| >0 if lk(t)||>Pvi, (2.7) 

Thus a configuration space trajectory t i— > of energy leaving the ball 
ZZ(E) C M d of radius R v i T (E) (the interaction zone) cannot reenter ZZ(E) in 
the future but goes to spatial infinity. Namely assume that (g(0),p(0)) > 0. By 
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so that 

||?(t)H 2 >lko|| 2 + l^ 2 (t>0). (2.8) 

Thus after having shown existence of a flow $ G C 1 (M x P, P) we can use 
results, derived in |DG] for scattering by non-singular potentials, in the singular 
case, too. In particular, we have 

limsup ||g(i, x)\\ — oo if and only if lim \\q(t, x)\\ — oo, 

and similarly for t — > — oo. 
For E > Hill's region 

n E :={qEM\ V(q) < E} 

is non-empty, but need not be connected (since there may be potential pits). 

By the assumption limiuii^oo V (q) = 0, IZe contains the neighbourhoods of 
infinity of the form {q e M \ \\q\\ > R} for R > large. These are connected 
if and only if d > 2. So for d > 2 there is precisely one unbounded connected 
component 1Z% of TZe, and the same is true for the energy shell projecting 
to Hill's region. We denote the unbounded connected component of by Eg. 

Definition 2.3 • We call E > a non-trapping energy if no orbit in is 
trapped, that is, for no initial condition x E T, E 

lim \\q(t, x)\\ = oo but limsup \\q(t, x)\\ < oo. (2-9) 

• The set of non-trapping energies E 6 (0, oo) is denoted by NT . 

Remarks 2.4 1. Unlike in Def. O in Def. 2.1.3 of [DS], E is called trap- 
ping if there exist orbits in bounded at least in the future. 
That definition has some advantages in the context of semiclassical quan- 
tum mechanics, but it would unnecessarily narrow the scope of our results. 

2. As $ is reversible, Def. 12.31 does not change under a sign change in (12.91 ). 

3. Trivially trapped and scattering orbits only occur in whereas the orbits 
in T,e \ are bounded, but there may be bounded orbits in Y, E as well. 

4. As shown in Prop. 1 of |Knj . in S 1 ^ existence of trapped orbits and existence 
of bounded orbits are equivalent properties. 

In particular for E £ NT there is no rest point in S^. But this implies 
that E is a regular value of H on S^, so that T, E is a smooth manifold. 
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5. The set NT C (0, oo) of non-trapping energies is open (see the proof of 
Prop. 2.4.1 of [DG]). 

6. As an example of physical relevance, for Coulombic potentials of the form 

n „ 

k=i Sfe " 

in the repelling case (Z fc < 0) there exists an interval (0, E ) C NT (see 
Sect. 5 of [OK] ). 

For n > 2, independent of the signs of the charges Z k , for s 1; . . . , s n G M 3 
in general position there exists an interval (E'th, oo) of trapping energies, 
where the dynamics of the bounded orbits is homeomorphic to the one on 
the suspended flow for the shift space (II. lft (see |Kn2] , Thm. 12.8). 

For E G NT the asymptotic directions 



and impact parameters 

<7.f : - < ?f(a:) := Jjm^ (<?(*, x) - (q(t, x),p ± (x)) p ± (x)^ 

are continuous ^'-invariant functions (see [DG], Thms. 2.5.2 and 2.7.2). 

By its definition, the impact parameter is orthogonal to the asymptotic di- 
rection, and for non-trapping energies E G NT we obtain homeomorphisms 

A% : - iV := T*S d ~ l , [z] ~ {qf{x),p±{x)). 

between the space of unbounded orbits and N . For E G NT the scattering map 

S E = (Qe, Pe) := A% o (A^)- 1 : N - N (2.10) 

is a homeomorphism of the symplectic manifold (N,un) and in fact a symplec- 

tomorphism, as follows from [DG], Thm 2.7.11. 

In particular for each initial direction 9 G S^ 1 the restriction 

E,e ■ i-e ^ o 

of the final direction map Pe ■= itsd-i o Se ■ N — >• is continuous. 
Lemma 2.5 For all dimensions d > 2, energies E > and directions 9 G S"^ 1 

lim P E Al±) = d - 

\wx\\-*oo 
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Proof. By continuity of (q2,P~) : £^ — > iV = T*5* d_1 and compactness of the 
interaction zone TZ(E) there is an R > with the property that for \\q±\\ > R 
the orbit (A E )~ 1 (q±,9) C T. E does not intersect the compact in Y> E lying over 
1Z(E). 

Thus there is exactly one point x(q±) = (po, 9o) £ Eg on that orbit whose 
configuration space projection q has minimal norm. By increasing the above R, 
that minimal distance diverges. As (po,qo) = 0, similar to ( 12.8ft we have the 
estimate 

||g(t)f> \\q4 2 + \Et 2 {t G E) 

for the w//?o/e trajectory. Integrating the force field, that is, the negative acceler- 
ation, along the trajectory we get uniformly on T, E , by using ( 12.2ft 

/ \\F(q(t))\\dt^0 for ||g || -> oo. 

Thus in this limit the change of velocity and of direction go to zero. □ 
By one-point compactification (t^S^ 1 U {oo}) = S^ 1 of that (d - 1)- 
dimensional vector space we may thus extend Pe s uniquely to a map 

P E ,e ■ S*- 1 -> S^ 1 (E g A/T, 6 g S^ 1 ). (2.11) 

which is jointly continuous in its argument and parameters. The choice of an 
orientation on the sphere fixes an orientation of the cotangent space TqS^ 1 , 
too, and we denote by 

deg(E) := deg(P^) 

the topological degree of this map. 

In general the degree of a map / G C 1 (S' d ~ 1 , S^ 1 ) is given by 

deg(/) = sign det (£>/(£)), 

xef-Hv) 

evaluated at an arbitrary regular value y of /. Then this definition is uniquely 
extended to C(S d ~\ S^ 1 ) (see, e.g., Hirsch [Hi], Sect. 5.1). 

In our case the degree is independent of the choice of orientation on S' d_1 . 
By joint continuity of Pe,8 m its argument and parameters it is also independent 
of the choice of initial direction 9. So the non-trapping degree 

deg :NT^Z 

is well-defined and locally constant on the (open) set of non-trapping energies. 

In |Kn] the degree was calculated for centrally symmetric (V(q) = V(\\q\\)) 
potentials, with the following results for regular values E of V: 



7 



For non-singular V 



deg(E) 



1 ,Ee(0, V max ) 
, E e (F max , oo) 



Here &R,% is homeomorphic to a (d — l)-sphere if E E (0, V max ) and 
dft| = for EE (V max ,oo). 

For singular V of the form V(q) = —Z\q\~ 2n ^ n+1 ^ all energies E > are 
non-trapping and 



deg(£) 

This is illustrated in Figure [TJ 



— n , d even 

;(!-(-!)") ,rfodd 



(2.12) 




Figure 1: 2D scattering with degrees 1 (left), (center) and -1 (right) 
In the present article we will show that 

• also in the case of general non-singular potentials V only the two cases 



and dlZ% = S are compatible with E being a non-trapping 
energy (Thm. loTi) . and that in this case deg(^) coincides with the above 
values (Thm. 13. ID . 



• Likewise, it will be shown in Thm. 14.31 that the degree formula (12. 12H 
remains true for arbitrary smooth perturbations of the singular potential V 
and energies E e N'T . 

We set the stage by parametrizing the sets 



J E,b 
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of phase space points in T, E with initial direction 9. By Theorem 2.7.2 of [DGj, 
given initial values x ,x G L Ee , the limit 

lim (q(t,x) - q(t,x )) 

exists, and there exists a unique time T Xo (x) G K such that 

lim (q(t + T X0 (x),x) - q(t,x ),p~~(x )) = 0, 

t— » — OO 

thus asymptotically synchronizing the two trajectories. 

For all gj_ G TqS^ 1 there is a unique phase space point X = X(g ± ,6*) in 
L E with 

lim X) - x )) = q±. 
This gives a parametrization 

/ 9 :RxTXUE» , (UJk^oI^^) (2.13) 

of the embedded C 1 -submanifold L E0 C T> E by a d-dimensional vector space. 
L Ee approximates for t — > — oo the affine Lagrange space 

{(p,q) e T * M I p= v 7 ^ 

in the C 1 -sense 

lim p(t,X(q ± ,9)) = V2E 6 , lim D (tiqx) p(t,X(q ± ,6)) = 

(see [DG], Thm. 2.7.1.). Since L Ee is also invariant under the symplectomor- 
phisms $*, it is thus a Lagrange manifold, too. 

3 The Case of Regular Potentials 

In this section we consider potentials V G C 2 (lR d ,lR) satisfying the long range 
estimate ( 12.20 . and energies E > 0. Then V~ 1 (E) C dlZ E for the unbounded 
component TZ E of Hill's region. 

Theorem 3.1 For non-trapping energies E G NT of V the following holds 
true: 

1. ifdTZ u E = and d > 2, then deg{E) = 0. 

2. if&R E = S d -\ then deg(E) = 1. 
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Remark 3.2 Note that only these two cases can arise for non-trapping energies, 
see Thm. 16.11 



Proof. • In case 1) we have \\p\ = y/2(E - V(q)) > for all (p,q) G T> E and 
thus use the continuous map (with x A := x/||x||) 

p-.^E^S*- 1 , (p,q)^p\ (3.1) 

and (similar to (12.131 )) the reparametrized flow on with initial direction 9 

Ie-.TgS^X (-1,1) , (q ± ,t)^^(tanh- 1 (t),X(q ±1 6)). 

We uniquely extend their composition 

h : TgS^ 1 x (-1, 1) -> S^ 1 , ~h = poI e 

to amap/i: S^ 1 x [-1, 1] -> S d -\ with compactification (t^S 4 - 1 U {oo}) = 

h{oa,t):=6 , h(x,-l):=9 and /i(a;, 1) := P E ,e(x). (3.2) 

By Lemma l2~5l /i is continuous. So the restriction maps h\ S d-i x s i \ (i = —1,1) 
are homotopic. Thm. 1.6 of [Hj], Sect. 5.1, gives the middle equality in 

deg(E) = deg (h\ S d-i x{1} ) = deg (h\ sd -i x{ _ 1} ) = 0, 

whereas the right equation follows from constancy of h\ S d-i x rn- 
• In the second case (371% = S^ 1 ) the statement follows for d = 1 trivially 
from the definition of the degree of a map S° — > S°. So we assume d > 2. 
Like in (13.20 . we construct a homotopy h which interpolates between Pe,6 and 
an explicit map on the sphere whose degree we can determine. 

Now (13.11) is not defined everywhere, and we use the bounded smooth map 

p : S s -> R d , (p, q) h-- 



2E 

instead of p. Note that \\p(p, q) —p(p, q)\\ — > uniformly in (p, q) as \\q\\ — > oo. 

Without loss of generality, we assume that G IR d \ T^. Then we have a 
smooth (C 2 ) map 

d-1 i \ , , „A _ Q 



q-.Zl^S*- 1 , (p,q)^q 

and unconditionally lim^oo q o <S>(t, x) = p + (x), with asymptotic direction p + . 
Unlike in case 1), we interpolate between q and p, using as parameter the con- 
tinuous distance function 

Dist r : -> [0,1) , (p,g) h+ tanh (inf {||g - cB\\/r \ c < 0}) (r > 0) 
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from the incoming axis defined by the initial direction 9 e S^ -1 , with a suitable 
parameter r. 

More precisely, with the impact parameter w.r.t. the initial direction 9 

q ± ,e:^ U E-T e *S d - 1 , (p, q) -> q - (q,9) 9, 



(x) := xj a/||^|| 2 + 1 and i_ := max(— t, 0) we define the map 
h : T* e S d ~ l x (-1, 1) -> S^ 1 

(Dist r • p + (1 - Dist r ) -(q + t- (q±,e)) A ) o I e (x, t) 



h(x, t) 



(Dislv • P + (1 - Distr) '(q + t- (q±,e)) A ) ° f) 



• We begin by showing that this is well-defined for r large. 

First, q + 1_ 7^ since in general || (x) || < 1 so that || < 1. 

Thus the numerator and denominator of h are well-defined. 

The denominator could only vanish on the hypersurface defined by Dist r T 1 (l/2). 

Projected to configuration space, this consists of the union of 
the hemisphere {q G R d | \\q\\ = rtanh _1 (|) , (q,9) > 0} 
and the cylinder {q e M. d | \\q - (q,9)9\\ = rtamV^i) , (q,9) < 0}. 

- For (p,q) = h{x,t) projecting to the hemisphere we have (9,q) > 0, 
{9,q±,e) = 0, and for all e > (9,p) > —e if r(e) is large. The last 
statement follows since for orbits in 

- not intersecting a ball of radius R > R vir in configuration space the 
change of maximal momentum is uniformly going to zero as R — > 00, 

- whereas for the orbits intersecting in configuration space that inter- 
action zone and then the hemisphere of radius r > R at (p,q) G Eg 
the difference \\p — q\\ uniformly goes to zero as r — > 00. 

- On the cylinder we have 

- a similar inequality for the outgoing parts of those orbits which have 
intersected the interaction zone. 

- For the incoming parts 

{$(t,x) I q((-oo,t],x)nlZ(E) =(&} 



of the orbits 



(9,p)>l-e and (9, (q + q ± , e ) A ) = ..^ „ > -1/2, 

\\q + q±,e\\ 

since for radius r > 1 the denominator \\q + q±,$\\ > 2. 
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A 



So in both parts of the hypersurface the denominator of h does not vanish. 
• We now consider the limit behaviour of h and define an extension 

h : S*- 1 x [-1,1] -> S**- 1 

of h using these limits. 

- h(x, 1) := linis^! s) = P E ^_(x). Namely, in the large time limit £_ 
vanishes, and lim^ip o J e (x, s) = lim s _^i g o Iq(x, s). 

- h(oo, s) := limii^n^oo h(x, s) = 9 (s G (0, 1)), 
since then lini|| :B ||_ ) . 00 Dist r o Ig(x, s) = 1. 

- The limit h(x, —1) of early times is given by 

lim h(x, s) = (tanh (^Pj Q + ^1 - tanh (^Pj J (x - 9) A 

This is a continuous map T^^ -1 — > 5' d_1 . After compactification to a 
continuous map S^ 1 — > S^ 1 (setting h(oo, —1) := 9) it is of degree one, 
as follows from linearization at the unique preimage of —9. 

So similar to 1), by one-point compactification of the cotangent spaces we can 
uniquely extend h to a continuous map 

h : S*- 1 x [-1, 1] - S*- 1 , with h\ sd - lx{1} = P E> g. 

The two restriction maps h\ S d-i x u\ are homotopic so that 

deg{E) = deg (h\ S i-i x{1} ) = deg (h\ sd -i x{ _ 1} ) = 1. 

This shows the validity of the second claim. □ 

Remarks 3.3 1. For the class of regular potentials meeting the inequality 

(q,VV(q))<0 (qeR d ), (3.3) 

J\fT D {E > | E regular value of V}, as can be seen by comparison 
with the virial identity (12.51 ). 

2. We may also consider scattering by an obstacle B C M d diffeomorphic 
to a d-dimensional ball. Then OB = S^ 1 and we have the Gauss map 
n : dB — > S^ 1 . We use the cotangent bundle T*(dB) of the boundary 
of the obstacle to describe the reflection data. Using the euclidean metric 
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on configuration space, T*(dB) is considered as the 2(d — l)-dimensional 
submanifold of T*M. d annihilating n. 

Without loss of generality we fix the value of the Hamiltonian 

H:T*M^R , H( Pl q) = ±\\p\\ 2 

on M := R d \ B to be \ so that \\p\\ = 1. This implies that the tangential 
component p — (p, n(g)) n(g) of the momentum p is contained in the unit 
disk of T*(dB). Then scattering means to invert the normal component 
and to leave the tangential component invariant. 

As that normal component vanishes for solutions tangential to B, the 
scattering map is still continuous (though not continuously differentiable). 
So if we assume that the obstacle is non-trapping, then we get index 1, by 
the same argument as in the above theorem I3~T1 As an example, the non- 
trapping condition is met if B is star-shaped, since then every reflection at 
a point q G dB increases the value (q(t),p(t)), in analogy to (13.31) . 

4 The Case of Singular Potentials 

We now treat the case of a singular potentials of the form 

V(q) = - r, Z lla +W(q) (a 6 (0,2)), (4.1) 

\\q ~ s \\ 

introduced in Sect. |5] Due to the singularity at the origin the Hamiltonian flow 
in the phase space T*M is incomplete. However, for certain values of a this flow 
can be completed by phase space extension. Then the regularization is essentially 
unique. 

For different regularization schemes of the representative Kepler problem 
(with or without time change and change of phase space dimension) consult 
Chapter 11.3.4 of Cushman and Bates [CB], and Chapter 5 of Cordani [C"o| . 

Proposition 4.1 For d > 2 the following statements are equivalent: 

1. a = 2n/(n + 1) for n e N. 

2. The phase space T*M can be extended to a 2d-dimensional symplectic 
manifold (P, uj), with a C 1 flow $ : R x P — > P, extending the incomplete 
Hamiltonian flow generated by ( 12. 4|) . 

Moreover in this case P is a union P = T*M U(Ix S^ 1 ) and can be given the 
structure of a 2d-dimensional symplectic manifold (P,u>), extending (T*M,u> ). 
H : T*M — > R then extends to a function in C 2 (P,R), also denoted by H, 
having the same regular points, and its Hamiltonian flow is $ : R x P — > P. 
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Proof. • We denote the maximally extended (incomplete) Hamiltonian flow on 

T*M by $ : D -> T*M, with open domain D C R x T*M on extended phase 
space. As follows from general theory of o.d.e., D is of the form 

D = {(t,x) GRx T*M \tete (T~(x),T + (x))} 

with escape times T~ : T*M — > [— oo,0) := {— 00} U (— 00, 0) upper semicon- 
tinuous and T + : T*M — > (0, 00] lower semicontinuous. 

By reversibility of the flow we consider only T + . Like in Thm. 3.1 of [MG] 
we conclude that for T + (x) < 00 we have a collision at time T + (x), that is 

lim q(t, x) = s. 

t/T+(x) 

Without loss of generality we assume s = 0. Similarly we assume that Z = 1, 
using a rescaling. 

• We first assume that W = in (14.11) . Then due to the centrally symmetric 
form of V every solution curve t 1— > q(t,x ) in M with initial conditions x = 
(po,<?o) ^ T*M lies in the plane (or line) spanned by p and q . 

• So we can assume d = 2 for the moment. The angular momentum 

L : T*M — > R , L{p,q) =q lP2 -q 2Pl 

is conserved by the maximally extended Hamiltonian flow $ : D — ► T*M. For a 
trajectory with energy E and value / of L we calculate the total deflection angle 
A<p(E,l), as seen from s. 

Considering for a moment an arbitrary centrally symmetric potential V(q) = 
V(\\q\\) and for I ^ its effective potential Vi (with VJ(r) := V(r) + J^), there 
may or may not be a largest r > with Vz(r) = E, then called the pericentral 
radius r min . In this case we have (see Chapter 2.8 of Arnold [Ar]) 

f°° Lb f°° l/r 2 

Acp(E,l)=2 -hdr = 2 ' =dr. (4.2) 



r 



^(E-V+r)) 



Setting V(r) := —r a with a E (0,2), we see that r min is well-defined and 
non-zero for / 7^ 0. Substituting v := OUAv^L we obtain 

/""max J„. 

A<p(E,l) = sign(l)2 . (4.3) 

1— a 



2£|Z| 1 -«/ 2 2^ +w Q - 

with 2E\l\ T =^2^ + f^, ax -f^ ax = 0. Since a < 2, in the collision limit / -> 
the first term in the square root vanishes, and 

f 1 dv 9-tt 
Alo ± : = lim A(p(E,l) = ±2 / = ±_f!!_, (4.4) 
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which equals ±(n + l)n if a = 2n/ (n + 1). 

So precisely for those exponents a G (0,2) that appear in our first assertion we 
have A(/9 + = Ay?~ (mod27r). This shows the implication 2) =^> 1). 
Moreover for a = 2n/(n + 1) formula (14.3ft equals 



A<p(E,l) = sign(Z)2 



dv 



V / 2 J E/ 2n 2( 1 " n )/ 2 + W 2n/("+l) 



2 



which with (14.4ft extends to a 5 1 -valued function of I and _E which is smooth 
even at / = 0. 

• In order to prove the implication 1) =>- 2), we now assume a = 2n/(n + 1) 
with n E N. Then we can continuously regularize the collision orbits with I = 
after collision at time to, by setting 

(p(to+At), g(t +At)) := ((-l)>(t -At), (-l)" +1 g(t -At)) (At > 0). 

Still that trajectory is undefined for time t , since q(t ) M. 

• P can be made a 2d-dimensional manifold and $ a smooth Hamiltonian flow, 
by using adapted coordinates in a suitable phase space neighbourhood 

\\q-s\\ < e, ||pf > 1 



U £ := { (p,q) e T*M 



\\q-s\\> 

with c a := 2 -1p G (a, 2). For small £ within t> £ 

|<?-s,P> = \\p\\ 2 -j^^-(q-s,VW(q)} 

> Ca 2 a \\ q -st - {q - S > VW{q))> °- (4 - 5) 
So within U £ the flow is transversal to the pericentric hypersurface 

S :={(p,q)eT*M\ {q-s,p) = 0}. (4.6) 

As Jjzll? — s|| 2 = 2^ (q — s,p), this inequality also shows that the point of the 
orbit on So is indeed pericentric. Every collision orbit enters U £ , since 

- ° aZ -2^4^ + 2(E-W(q))^oc 



\\q-s\\<* \\q 
as q approaches s. 

• In the present case W = we use the following coordinates on U £ 
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H\ tjE G C°°(U £ ,R). The value E := H(x) of the Hamiltonian function 
at x is conserved by the flow. 

The time T : U £ —> R needed to arrive at the pericentre respectively at 
s. As the Hamiltonian function is smooth, we have $ G C°°(D,T*M). 
Furthermore by (14.51) the flow $ is transversal to the smooth pericentric 
hypersurface So, so that T is smooth for all points x := (p,q) G U £ on 
non-collision orbits. 

Moreover, T is explicitly given by the integral over inverse radial velocity: 

r\\4 i 
T(x) = sign((p,g)) / - ^= dr, 

J r ^„ ./2(E-Vi(r)) 



sign((p,g))2" 



p(n+i) dv 



t> 2 v / 2£'Z 2n 2( 1 ~ n )/ 2 + v 2n /( n + 1 ) — v 2 

2 (n + l)/2|| q || 



with E := H(x), I := L(x), and similar to A^, T G C°°(U £ , R). 

As in [MG] , we now discern the cases of even resp. odd n G N. 
In each case we define a map F G C(U £ , S^ 1 ) in a way so that for initial 
conditions x = (PO)9o) G U £ the direction F(xo) lies in the two-plane (or 
line) spanned by p and q — s. 

For odd n and non-zero angular momentum of xq we define F{xq) G S^ 1 
as the direction of the pericentre of the orbit through x . For t := T(x ) 
(that is, $'°(a;o) G S ) this equals 

q(t ,x ) - s 



lk(*o,^o) - s|| ' 



Dividing the expression (14.40 for the limiting deflection angle of collision 
orbits by 2, we note that for zero angular momentum of xq and for a = 

2n/(n + 1), n odd we get F(x ) = (-l) (nH 



11/2 go- 



ll<?0-s|| ' 

For even n and non-zero angular momentum of xq we define F{xq) as 
the normalized velocity F(xq) := n|^^j] at the pericentre of the orbit 
through x . By (14.61) this is perpendicular to the vector q(t ,x ) — s. 
So using formula (14.41) for the limiting deflection angle of collision orbits 
we note that for zero angular momentum of x and for a = 2n/(n + 1), 
n even we get F(x ) = (-l) n/2 ^- 

The conserved (non-zero) value I of the 'angular momentum vector at 
the pericentre' L(x) := ||g(t 0) ^) — s\\p(t ,x). (l,<p) is a point in the 
symplectic manifold T*S d ~ 1 . 
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The collision orbits correspond to the points with I = 0, but T(p,q) ^ 0. The 
cylinder K. x S^ 1 in ( 14.7 ft is then identified with the set of missing phase space 
points, characterized by (l,t) = 0. 

In the above coordinates the flow is affine in the variable T (T o $* = T + £)), 
the other variables being constants of motion. So $* can be uniquely extended 
to the cylinder, and the resulting flow on P is smooth and complete. 

That collision orbit can thus be parametrized by its energy E 6 K. and, say, 
initial direction 9 E S^ 1 . So by setting 

P:=T*M0(Rx S*" 1 ), (4.7) 

we may thus regularize the motion on this new phase space and obtain a complete, 
smooth flow extending $ 

$ : R x P -> P. 

• If the smooth potential 14 7 in (14.1ft is non-zero, the above quantities H,L and 
F are not conserved. However, they can be used to define conserved quantities, 
namely their values at the unique pericentre of the near-collision orbit. 

See [KK], Prop. 2.3 and |Kn2] Thm. 5.1 for details of the (somewhat tech- 
nical) construction in the representative case n — 1 of the Kepler potential. 

• We now extend the natural symplectic form on T*M to P, defining it by 
($*)*cl> on P \ T*M. More precisely, by (03]) for any x E P \ T*M there is an 
open neighbourhood U cPofx and t > such that V := $(t,U) C T*M. 
The restriction * := ¥\ v is u; -symplectic on V := n T*M). We 
uniquely extend cu to P by setting cu\ u := ^*cu (Concrete expressions of to in 
terms of local coordinates can be found in [KK], Prop. 2.3 and [Kn2] Thm. 5.1). 

• That the Hamiltonian function extends to a function H E C 2 (P, M) having no 
singular points on P \ T*M and generating $, is immediate from the forego- 
ing construction, since H is one of the coordinates used in the definition of P. □ 

In the regularizable case for E > V mSiX the energy surface Y> E is a (d — 1)- 
sphere bundle 

7T E : S £ - R d (4.8) 

over configuration space. As the base M. d is contractible, this bundle is trivial. 
The same statement applies for all directions 9 to the induced bundles 

& ■= (he o IeT^E :t E ^Rx T* e S d ~ l (9 E S^ 1 ) 

over the (parametrized) Lagrange manifolds L Ee , see (12.131) . By definition of 
induced bundles (compare with [Hi], Sect. 4.2) the total space of £ 9 equals 

t E := i;X E := {(x,y) Glx T^S^ 1 x Y, E \ I e (x) = ir E (y)}. 
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However, if we consider the local degree e(S) of the section 



S:M^E E , , (y/2(E-V(q)) 6,, 

see Bott and Tu |BT] , §11, then this is non-trivial. By definition this is the 
degree of the composed map 

S^ 1 dB r H E \ Br ^B r x S*- 1 ^ S*' 1 , (4.9) 

for a ball B r := {q e M. d \ \\q — s\\ < r} of arbitrary radius r > 0. 

Remark 4.2 In the case of sphere bundles over a compact base manifold M, 
the Euler number of the bundle is the sum of local degrees at finitely many base 
points, see [BTJ, Thm. 11.16. 

In this context the Euler number vanishes for d = dim(M) odd, if the oriented 
(d — l)-sphere bundle is the restriction of a vector bundle of rank d over M. 

As Formula (14 . 10D below indicates, for our bundle tie '■ S B — > M. d this is not 
the case if n is odd. 

Theorem 4.3 For a non-trapping energy E E NT of V the following holds 
true. IfV is of the form ( TO) with a = 2n/(n + 1) for neN, then for d > 2 



i — n , d even 
deg(E) = e(S) = \ d Qdd (4.10) 



Proof. • We start by calculating the local degree e(S). Without loss of generality 
we assume that the singularity is located at s = 0, with a = 2n/(n + 1). 
• We consider first the case of a potential (12. lft with W = 0. Then, as there 
are no (semi)-bounded orbits of positive energy, NT = M + , and we can use the 
formula 

Ay? = \imA(p(E,l) = nir (4.11) 

for the limit of the total deflection angle, derived in (14.4ft . 
We define a section of the bundle 7Te : Eg — > M d by 



T : 



T (q)-=l (y/W-V(q))F(0,q),q) ,q^s 
I (E,F(B,q)) ,g = s 



where, similar to the proof of Prop. 14.11 F : T. E — > S^ 1 maps x to the unique 
pericentral direction of the orbit through x. This section is continuous, the 
apparent discontinuity of (14.12ft at q = s being owed to the use of the local 
cylinder coordinates (E,(p) elx S d ~ x . 
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• Evaluating the degree of (14.91) in the limit r — ► oo, we can use (14.111) to obtain 
the second equality in (I4.10p . 



In the case of d = 2 dimensions the outgoing angle 

PeAV-l) = 6 ~ A V(E, V2Eq ± ) (q± > 0) 
is continuous decreasing in q±. So in this case it follows from (14.40 that 

f°° d 

/ -t-VeMa-) d Q± = "2Ay = -2™. (4.13) 

J -00 dq± ' 

This is twice the change in direction from 9 to F(9,q), since by symmetry the 
change in direction before and after the time of pericentre are equal. The section 
T : M. d — > Tje based on F trivializes the circle bundle, and on dB r the difference 
between the sections S on T is given by g 1 — > F(9, q) — 9. 

On the other hand half of (14.131) is the contribution of the part {q E dB r \ 
(q,9) < 0} (on the left hand side of (14.9ft ) to e(S), since it corresponds to 
the incoming parts of the orbits. By symmetry the outgoing parts of the orbits, 
corresponding to {q E dB r \ (q,9) > 0}, give the same contribution. 

Together this proves 

deg{E) = -n (E > 0). 

• For d > 2 we consider a family of trajectories with fixed E and 9, whose impact 
parameter q± varies on a one-dimensional subspace L C T(jS d ~ 1 . 

9 and this subspace span a 2-plane in W l , and 9 + lies in that plane. To avoid 
degeneracies we choose a 9 + which is linear independent from 9. Then there are 
exactly n impact parameters q]_, . . . , q™ E L with P E ,e(q±) = 9 + . 

[n/2] of them have a scalar product (q]_,9 + ) > 0, and (q]_,9 + ) < for the 
rest. For the first group the restriction of the linearization of the final angle 
map to the subspace {v E TqS^ 1 \ v J_ L} gives a positive sub-determinant, 
whereas for the second group the sign equals (— l) d ~ 2 . So 

deg(E) = - ([n/2] + (-l) d " 2 (n - [n/2])) , 

proving the second equality in (14. 10ft . For the case W = this also proves the 
first equality in (l4~IUi) . using ( 1Z121 . 



Now we turn to the case of non-vanishing W. The local degree e(S) is 
independent of the radius r > in ( 14.9ft . Evaluating e(S) in the limit r — > 0, 
we see that by smoothness of W it coincides with e(S'), calculated above for the 
case W = 0. This proves the second equality in (14.10ft for arbitrary W. 

The map (n E ,p + ) : — > M. d x is another trivialization of the bundle, 
and p + , evaluated over the sphere dB r x {9}, has the degree deg(E). This 
shows in the general case that e(S) = deg(E). □ 
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5 Projection of the Lagrange Manifold 



This section applies to regular as well as to singular potentials. 

As noted in Sect. [5J for E £ AfT and 9 G S"* -1 , the image L Ee of the 
embedding 

is a Lagrange manifold in phase space P. We now consider the projection n : 
P — ► R ' and the composition map 

II S := 7T o J fl : R x Tg*^ 1 -> 7^, 

mapping this Lagrange manifold to configuration space. This is a C 1 -map be- 
tween (i-dimensional <9-manifolds. Moreover, it is proper, that is, compacts have 
compact preimages. 

We orient the vector space R x T^S^ 1 so that for all x G W 1 ' 1 

det(DU E (t,x)) > for t<0 

(then, in fact, lim^-oo det(DU E (x, t)) = y/2E). 
Then for every regular value q of Tl E , we set 

deg g (n £ ):= sign{det{DU E {y))). 

yen- 1 ^) 

This is well-defined. By properness of II^ and connectedness of 1Z E , the value 
does not depend on q. Thus we obtain an integer 

deg(U E ) G Z. 



Theorem 5.1 deg(E) = 1 — deg(n E ). 

Proof. We evaluate deg (IIe) for q = H E (x, t) in the limit t — > — oo. Then, 
as remarked above, the point (x,t) G U E 1 (q) contributes a one to deg(IlE). 
The other contributions come from points (x', s) with s — > oo as t — > — oo and 
p + (X(x', s)) — > — 0. So they correspond to the terms contributing to 

deg(P E , e ) = 

sign det (DP E ,e{q±)), 

but the signs in the sum for deg(II^) are reversed, since the initial direction 9 is 
reversed for these orbits. □ 
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6 A Topological Criterion for Trapping 



This section applies to regular potentials. The following topological criterion 
for trapping generalizes a low-dimensional (d < 3) result of [Kn] to arbitrary 
dimensions d. 

Theorem 6.1 If E e NT for a long-range potential V : R d -> R, d G N, then 
the boundary dlZ E of Hill's region is empty or homeomorphic to S^ 1 . 

Proof. • If E G NT, then the relative homotopy groups of Hill's region w.r.t. 
its boundary are trivial, 

ir k (ni,ani) = {e} {ke{o,...,d}). (6.1) 

This was shown in Thm. 3.2 of [Kn] for short range smooth potentials, but the 
argument only involved the dynamics within the interaction zone XZ(E) and thus 
generalizes to our class of long range potentials. 

Our aim is to invoke the /i-cobordism theorem in order to show that (16. ID 
implies d7Z E = or d7Z E = 5 d_1 , = denoting existence of a homeomorphism. 

• We assume that &R U E ^ and have to show d1Z u E = S"* -1 . First of all, for 
E G NT, dlZ E is a closed (d— l)-dimensional manifold, embedded in M d : 

1. E is a regular value of V\ n u for otherwise there would exist a equilibrium 
point of the flow, restricted to T, E . By Remark I2~41 4 this would contradict 
the assumption E G NT. As d1Z E is a component of V" 1 {E), it is a 
(d— l)-submanifold, and being itself a boundary, without boundary. 

2. Since by assumption limngn^oo V(q) = but E > and V\ &R u = E, the 
boundary dlZ E of Hill's region is compact, and thus closed as a manifold. 
We choose R > so that it is contained in the interior of the ball -B^(O). 

Moreover, we can assume dlZ E to be connected, for otherwise by curve shortening 
one could find a solution moving between two components of dlZ E (a so-called 
brake orbit, see Seifert [Se] and Gluck and Ziller |GZ] ). Again by Remark I2~4l 4. 
such a bounded orbit in T, E would be an obstruction to our assumption E G NT. 

• We want to apply the /i-cobordism theorem (see [Mi], Thm. 9.1) to the triad 
(W;V,072£), with 

W := ll E f] B d R (0) and V := dB d R (0). 

So the boundary of the manifold W equals <9W = VU dTZ"%. 

In fact the subset W of is homotopy equivalent to 72-^, as follows from 
the deformation retraction 

H : [0,l}xK E ^n E , H(t,x) := l m<R} (xyx+t m > R} (x)- (^(l-t)x+t„J 
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between the identity map 1^ and the map 



/ :TZ U E ^W , x i-> ^ii^ikh} (x) • x + 1 {M > R} (x) 




composed with the inclusion W — > 72.^. 

We can assume d > 4, since the cases d < 3 have been analyzed in [Kn], 
using low-dimensional methods. 

There are several assumptions to be checked in order to apply that theorem. 

• First we have to ascertain that W, V and d1Z E are simply connected. This is 
trivial for the (d — l)-sphere V. 

We claim that 7Ti(W, V) = {e} which then will imply = {e}, using 

the exact sequence 

{ e } = 7n(v) - m{W) - m{W, V). (6.2) 

A second implication will be that also &JZ E , the boundary of Hill's region, is 
simply connected. Here one takes the exact sequence 

{e} = ir 2 (K E ,dK E ) = ir 2 (W,dK u E ) - 7n{dTZ u E ) - tti(W) = {e}, (6.3) 

the first identity being (16.11) . the second using the deformation retraction H. 

• To show vanishing of 71"! (W, V) in (16. 2D . we construct for any two representa- 
tives hi : [—1, 1] — > W, hi(±l) E V of relative homotopy classes in 7Ti(W, V) a 
homotopy 

H : [0,1] x [-1,1] - W with H\ {iyx[ _ ltL] = hi (2 = 0,1). (6.4) 

This would be simple if we could assume that the representatives hi are solution 
curves of the Hamiltonian equations with energy E (with suitable reparametri- 
zations, mapping the time intervals onto [—1,1]). Then we could vary initial 
conditions from the one for h to the one for hi, and obtain a homotopy H. 

What we can do instead is the following. Without loss of generality we assume 
that the hi are smooth regular curves not meeting dlZ E . 

For the scaling constants c;(t) := ^/2(E — V(hi(t))/\\hi(t)\\ the initial con- 
ditions Xi(t) := (ci(t)hi(t), hi(t)) are in E^. The maps 

&:A->£jE , 9*{t,s) :=<S> s (xi{t)) 

with domains A := {(t,s) G [-1,1] x R \ s G [Tf(t), T t + (t)]} are C 1 - 
smooth, and the functions Tp < < are uniquely defined by the conditions 
gi(t,Tf{t)\ G V E for the disjoint hypersurfaces 

V±:={(p,q)Ei; E \\\q\\=R,±(p,q)>0} 
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over V. Namely, the flow line through Xi(t) G T, E intersects both surfaces, since 
E G NT. Furthermore by the virial inequality (12. 7 D both are intersected at 
most once by a flow line. By the same argument we get transversality of the 
intersection and thus C^-smoothness of Tf. 

As V E is connected, there is a homotopy ij) : [0, 1] x [—1, 1] — > V E with 

t) = $ (Tr(t), Xi (t)) (i G {0, 1}, t G [-, 1, 1]) . 

This gives rise to a continuous map p : D p ^ Y, E , p(z,t,s) := & s (tp(z,t)), on 

D p := {(z,t,s) G [0,1] x [-1,1] x [0,oo) \ s < S(z,t)}, 

with the continuous Poincare time S uniquely given by $(S(z, t),ip(z, i)) G V E . 
For a homeomorphism 

K : [0,1] x [-1,1] x [-1,1] -^D p , (z,t,s) ^ (z,t,k{z,t,s)) 

with k(z, t, -1) = 0, fc(z, t, +1) = 5(2, t) and t, t) = -T~{i, t) the config- 
uration space projection (see (14.81) ) 

F := ti e o p o K : [0, 1] x [-1, 1] x [-1, 1] -»• W 

has the property 

H(i,t) = hi{t) (i G {0, 1}, t G [—1, 1]) 

for H(z,t) := F(z,t,t). So is the desired homotopy from (16.40 . 

• As a second condition of the /i-cobordism theorem we have to ascertain that 
the relative homology 

H,(W,dn u E ) = 0. 

This follows from H*(1Z E , dlZ E ) = since W was shown to be homotopy equiv- 
alent to 1Z E . By the relative Hurewicz isomorphism theorem (|Sp|, Chapter 7.5) 
vanishing of that relative homology follows from (16.11) . 

• The last condition of the /i-cobordism theorem is that dim(W) > 6. In fact 
a version of the /i-cobordism theorem in the homeomorphic category works for 
d = dim(W) > 5, see Sect. 7.1 in the book |FQ| by Freedman and Quinn. 

In these cases its conclusion is that the following triads are homeomorphic: 

(W;V,dTZ u E ) = (V x [0, 1], V x {0}, V x {1}); 
in particular dlZ E is homeomorphic to a (d — l)-sphere. 

• The case not covered by the cobordism theorem is the one of dimension d = 4. 
Here the celebrated proof of the Poincare conjecture by Grigori Perelman allows 
us to conclude that &R% is homeomorphic to S 3 , since by (16. 3D dlZ E is a 
connected closed 3-manifold with fundamental group 7Ci(d1Z E ) = {e}. □ 
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